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ABSTRACT 


A  self-contained,  rapid,  Computer  Aided  Design  (CAD) 
program  for  a  desk  top  computer  (i.e.  HP  9845)  was  developed 
for  a  first  cut  approximation  for  the  design  of  a  super- 
cavitating  propeller  blade.  This  program  eliminated  the 
error-prone,  tedious  interpolation  of  empirical  data  graphs 
by  providing  approximations  and  curve  fitting  techniques  to 
augment  existing  formulae.  The  complex  Goldstein  function 
and  the  inexact  Prandtl  approximation  were  replaced  by  a 
more  simple  function,  the  Wilson  factor,  that  maintained 


the  confidence  level  of  the  manual  calculations. 
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I .  INTRODUCTION 


The  screw  propeller  was  first  developed  by  an  Englishman, 
Hooke,  in  the  17th  century,  but  received  little  use  until 
nearly  two  centuries  later.  In  the  late  1800's,  while  the 
British  steamer  TURBINIA  was  being  designed,  it  was  first 
observed  that  the  propeller  could  operate  in  a  cavitating 
mode  wherein  at  least  a  portion  of  the  propeller  operated  in 
an  ' air  cavity ' .  This  air  cavity  was  not  air  but  a  result 
of  critical  pressure  conditions  in  the  vicinity  of  the 
propeller  blades  [1] . 

Inasmuch  as  propeller  cavitation  usually  results  in 
marked  increases  in  required  rotational  speed,  propeller 
slip,  and  erosion  of  the  blade  surfaces,  as  well  as  decreases 
in  the  power  efficiency,  cavitation  was  to  be  avoided,  if  at 
all  possible.  One  precaution  was  to  insure  that  the  tips 
of  the  blades  were  always  well  submerged  to  as  not  to  have 
an  air  interface.  If  partially  submerged  then  ventilation 
may  result  yielding  similar  effects  as  cavitation. 

It  was  not  until  the  mid-twentieth  century,  while  search¬ 
ing  for  drive  systems  for  racing  boats,  that  it  was  discovered 
that  by  operating  in  extreme  cavitation  regions,  high  water¬ 
craft  velocities  could  be  obtained  with  but  minimal  deleteri¬ 
ous  effects.  In  fact,  not  only  could  the  propeller  be 
operated  in  a  cavitating  mode  fully  submerged,  but  it  could 
also  be  operated  with  a  portion  of  the  blades  out  of  the 


water  partially  submerged.  As  a  result,  these  systems  are 
referred  to  as  either  fully  submerged  or  partially  submerged 
propellers.  With  either  propeller  the  effect  is  essentially 
the  same  although  the  loads  may  differ. 

Whereas  the  design  techniques  for  conventional  propellers 
had  become  well  established  [1] ,  few  techniques  were  available 
for  the  supercavitating  blades  or  foils.  In  general,  the 
subcavitating  design  procedures  were  used  as  an  approximation, 
but  this  left  much  to  be  desired.  In  the  mid- 1950's,  Tulin 
[1]  proposed  two-dimensional  lift-drag  characteristics  for 
these  designs,  and  at  the  same  time  Morgan  and  Tachmindji 
applied  circulation,  or  lifting  line,  theory.  The  latter 
techniques  were  based  on  the  works  of  Prandtl,  Hemhold,  Gold¬ 
stein  and  Lerb,  utilizing  hydrodynamic  principles.  Now, 
nearly  thirty  years  later,  these  approaches  are  still  used  in 
the  design  of  supercavitating  propellers  {1] . 

The  most  rigorous  of  these  techniques  is  Lerb's  [2]  indue- 
tion  factor  method,  while  the  use  of  the  Goldstein  function, 
though  less  precise,  compares  quite  favorably  to  Lerb's  method. 
Any  of  the  methods  used  in  the  design  of  supercavitating 
propellers  have  been  labor  intensive  [2]. 

A.  STATE-OF-THE-ART 

The  formal  design  methods  for  high  speed  naval  ship  pro¬ 
pellers  were  established  by  Morgan  and  Tachmindji  of  the 
David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
(NSRDC) , [6] ,  using  the  Goldstein  function  to  correct  for 
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supercavitating  operations.  The  major  problem  in  these  tech¬ 
niques  was  attributed  to  blade  cavity  effects,  and  a  correction 
in  the  form  of  an  accountability  factor  is  used  to  offset  this 
effect.  This  correction  is  based  on  an  approximate  cavity 
thickness  in  the  lifting  surface  calculations. 

Another  approach  in  the  design  method  is  that  taken  by 
Hydronautics ,  Inc.,  [1]  that  utilizes  lifting-line  theory  and 
the  Goldstein  function  together  with  a  camber  correction  fac¬ 
tor  and  two-dimensional  foil  analysis  about  a  free  surface. 

A  comparison  of  the  two  approaches  by  Tulin  indicates 
an  overestimation  of  efficiency  and  thrust  by  the  NSRDC  method 
[1] .  Both  methods  rely  heavily  on  the  use  of  empirical  data 
and  require  frequent  interpolation  and  extrapolation  of  data. 
Most  of  the  more  recent  research  in  this  field  has  been  con¬ 
ducted  between  1955  and  1975,  and  the  confidence  level  in  the 
selection  of  a  supercavitating  propeller  without  full-scale 
testing  has  been  unfavorable. 

B.  FUNDAMENTALS  OF  CAVITATION 

The  phenomena  of  cavitation  is  the  formation  of  bubbles 
about  a  foil  or  blade.  This  can  be  explained  by  the  Bernoulli 
equation,  assuming  inviscid  compressible  behavior,  where  the 
sum  of  the  static  pressure,  p,  and  the  dynamic  pressure, 
ypV2,  is  a  constant. 

1  2 

P  +  |pV  -  pt  ■  constant  (1) 
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where : 


p  »  static  pressure  on  the  foil  downstream  side, 

V  =  velocity  of  the  flow, 

p  *  density, 

Pt  =  total  pressure  (total  head) . 

By  relating  the  constant  total  pressure  of  the  ambient 
conditions  (subscript  0)  to  the  constant  total  pressure  at 
any  point,  assuming  constant  density, 

p  +  JP  V2  -  p0  +  jpvJ  (2) 

By  combining  the  dynamic  pressure  terms,  it  can  be 
observed  that  the  pressure  at  any  point  may  be  expressed  as, 

p  *  pQ  +  Ap  (3) 

If  p  *  0,  so  that  pQ  »  -Ap,  the  fluid  changes  state  and 
cavitation  results.  In  reality,  the  static  pressure  does  not 
completely  reach  zero,  but  will  decrease  to  the  vapor  pressure 
of  the  fluid,  pv.  When  the  static  pressure  reaches  this 
value , 

P  -  Pv  •  P0  +  Ap  (4) 

the  fluid  'tears'  or  boils.  This  is  cavitation!  See  Figure 

1. 
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Figure  1. 


Picture  of  a  Propeller  Blade  Tip 
Cavitating  (NSRDC) 
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A  measure  of  the  amount  of  cavitation,  called  the  cavi- 
tat ion  index  or  the  cavitation  number,  is  determined  by 
dividing  the  pressure  change,  Ap,  by  the  dynamic  head  of 
the  flow. 

When  the  bubble  formation  that  commences  at  the  leading 
edge  extends  beyond  the  trailing  edge,  so  that  a  plume  appears 
to  form  downstream,  the  foil  is  said  to  be  supercavitating. 

C.  SUPERCAVITATING  PROPELLERS 

A  propeller  designed  for  optimum  non-cavitating  or  sub- 
cavitating  operation  cannot  perform  efficiently  in  a  cavitating 
environment,  and  if  it  is  desired  to  have  a  propeller  system 
that  can  utilize  the  high  speed  advantages  of  the  supercavi¬ 
tating  propeller,  a  special  design  of  the  foils  must  be 
accomplished  that  will  deliver  thrust  efficiently  in  a  quasi- 
gaseous  environment.  When  it  is  not  operating  in  a  supercavi¬ 
tating  environment,  there  is  no  real  advantage  to  this  type 
of  propeller  system,  as  shown  in  Figure  2. 

One  can  visualize  a  supercavitating  propeller  as  somewhat 
a  hybrid  operating  between  the  marine  and  the  air  environment, 
and  it  has  somewhat  the  appearance  of  a  stubby  propeller  of 
an  aircraft  rather  than  looking  like  a  conventional  marine 
screw. 

Current  design  techniques,  as  outlined  in  Chapter  II, 
call  for  extracting  values  from  empirical  curves  to  use  in 
the  calculations.  It  is  the  purpose  of  this  project  to  semi¬ 
automate  this  process,  as  outlined  by  references  1,  2,  and  3. 


II.  BASIC  DESIGN  METHODOLOGY 


The  criteria  for  use  of  the  supercavitating  propeller 
is  essentially  the  relationships  between  the  speed  of 
advance,  V  ,  the  rotational  velocity,  RPM,  and  the  cavita- 

OL 

tion  number,  a. 

As  in  any  design  problem,  the  first  requirement  is  the 
establishment  of  the  design  operating  point  for  which  param¬ 
eters  are  assumed,  known  or  desired.  Such  a  point  may  be  a 
speed  condition,  such  as  cruise,  hump,  maximum  or  patrol 
velocity,  or  it  may  be  a  power  condition.  There  are  several 
basic  parameters  that  are  unusually  well  known  or  can  be 
readily  assumed.  These  include  the  velocity  of  the  craft, 
resistance  (drag) ,  power  available,  speed  of  advance  (or 
wake  factor) ,  thrust  (or  thrust  deduction  factor) ,  number  of 
blades  and  submergence. 

As  a  first-cut  approximation  for  the  required  diameter 
of  the  propeller,  the  Burtner  equation  is  used  (1J, 

D  -  (50  x  (HP)°*2/RPM0,6]  ft  (5) 

The  thrust  is  non-dimens ionali zed  by  the  use  of  a  coeffi¬ 
cient  of  thrust,  Cj, 


(6) 
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where : 


T  «  thrust  (lbs) 
p  «  density  (slugs/ft2) 

D  ■  propeller  diameter  (ft) 

V  »  speed  of  advance  (knots) 

S  *  submergence 
A  *  disk  area  (ft2) 

The  speed  of  advance  is  also  non-dimens ionalized  to  an 
Advance  Coefficient,  J, 

J  -  101.34  x  V  /(RPMxD)  (7) 

A 

where  101.34  is  the  conversion  factor  in  terms  of  ft.  and  sec. 

The  above  factors  are  shown  in  pound-force,  pound-mass, 
and  feet  units,  but  the  entire  calculations  can  also  be 
accomplished  in  metric  units,  if  desired. 

Once  and  J  are  determined,  the  expected  efficiency  is 
obtained  from  empirical  data,  usually  plotted  as  J  versus 
/Cj^.  The  radial  distribution  of  pitch  is  determined  from 
Kramer  Thrust  Coefficient  Curves,  with  inputs  of  their  Advance 
Ratio,  A,  and  the  Ideal  Thrust  Coefficient,  C^,  where 

A  -  V  / (NxDxir)  (8) 

A 


and. 
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Qj.  a  1.07  Qj  (?) 

It  is  to  be  noted  that  only  approximate  values  are  used 
inasmuch  as  the  exact  value  of  CTi  is. 


Si  -  V(1-2eXi> 


(10) 


where : 

c  *  drag  to  lift  ratio 
A^  *  ideal  advance  ratio. 

At  this  point,  solutions  can  be  made  for  the  radial  pitch, 
(D  Xu  x\.)  ,  and  the  hydrodynamic  pitch  angle, 

Once  the  hydrodynamic  pitch  angle,  is  established, 

the  ideal  delivered  thrust  coefficient,  CTi,  may  be  deter¬ 
mined  by  using  the  more  rigorous  induction  factors  and 
numerical  integration  (Simpson's  Rule),  or  by  using  the 
Goldstein  function,  k,  [5] 


K 


z  r  n 

1  7T  U 

A 


(ID 


where : 

r  ■  circulation 

u  ■  axial  velocity. 

A 

The  delivered  ideal  thrust  coefficient,  C^,  is  compared 
to  the  required  value  of  CTi,  and  if  the  delivered  value  is 


f 

!< 

j« 


} 


l 
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less  than  the  required  value,  the  hydrodynamic  pitch  angle, 

0^,  is  corrected  and  the  problem  iterated  until  the  delivered 
value  of  is  equal  to  the  required  value.  If  the  iteration 
process  becomes  too  cumbersome  or  too  lenqthy,  interpolation 
may  be  used  to  validate  the  new  value  of  0^. 

At  this  point  the  designer  can  compute  the  geometric 
properties  of  the  blade,  first  by  computing  the  value  of 
Cl(£/D) ,  where. 


CL  -  L/(|pV2A)  (12) 

l  =  blade  width 

L  «  lift 


and  then  by  determining  the  blade  width  at  seventy  percent 
of  the  blade  span. 


CLlt/D)0.7xD 
°-7  '  *1 


(13) 


where  s 


C£  *  coefficient  of  lift  for  a  finite  foil. 
Reference  2  states  that  the  optimum  value  for  is  0.16 


and  that  value  is  used  in  these  computations 


A  recommended  distribution  of  section  chord  as  a  function 


of  the  chord  size  at  the  seventy  percent  span  [-2],  is  shown 
in  Table  I. 

TABLE  I 

RADIAL  DISTRIBUTION  OF  CHORD 


0.400 

1.088 

0.475 

1.086 

0.550 

1.073 

0.625 

1.046 

0.700 

1.000 

0.775 

0.903 

0.850 

0.763 

0.925 

0.565 

1.000 

0.000 

At  each  radial  station,  the  blade  shape  may  then  be 
determined,  as  shown  in  Figure  3. 


Figure  3.  Blade  Dimension  Sketch 
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From  the  hydrodynamics  of  lifting  surfaces  and  empirical 
data,  the  maximum  dimension,  (y /%)  ,  is  tl], 

ulctX 

(yA)„„  *  0.1553  CT  -  0.00462  a 

max  Xj 

where  the  angle  of  attack,  a,  is  a  function  of  the  lift 
coefficient. 

The  height-to-chord  ratio  is  then  corrected  by  using 
camber  correction  factors,  ^  and  k2,  are  functions  of  the 
expanded  area  ratio,  EAR,  which  is  the  ratio  of  the  expanded 
blade  area,  A^,  to  the  disk  area,  AQ, 

ear  -  VAo  *  ¥  J1  §  <**  (14) 

The  EAR  is  similar  to  the  BAR  blade  developed  area  ratio  with 
the  difference  being  that  the  EAR  is  with  respect  to  radial 
distribution. 

The  corrected  thickness  factor,  ly/Z)  .  is, 

c 

(y/*)c  «  kxk2  (yA)roax  (15) 

The  thickness  of  the  blade  section,  t,  is  determined  from, 
t/l  •  F'Cl  +  N»a  (16) 

where  F*  end  N '  ere  as  shown  in  Table  II  [1]. 
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TABLE  II 


THICKNESS  COEFFICIENTS  [1] 


d/Z 

F* 

N' 

0.0000 

0.0000 

0.00000 

0.000000 

0.0075 

0.01888 

-0.001326 

0.001254 

0.0125 

0.03244 

-0.002574 

0.001870 

0.05 

*  0.14189 

-0.016816 

0.005704 

0.10 

0.29150 

-0.039520 

0.010074 

0.20 

0.56636 

-0.078940 

0.017349 

0.30 

0.78458 

-0.097600 

0.022222 

0.40 

0.93187 

-0.092870 

0.024647 

0.50 

1.00000 

-0.071780 

0.025710 

0.60 

0.98340 

-0.037750 

0.025853 

0.70 

0.87797 

0.005940 

0.025605 

0.80 

0.68055 

0.056880 

0.024660 

0.90 

0.03886 

0.112090 

0.023047 

0.95 

0.20650 

0.139620 

0.022048 

1.00 

0.00000 

0.166950 

0.020937 

The  final  step  in  determining  the  geometric  parameters 
is  the  determination  of  the  corrected  pitch-to-diameter 
ratio,  P/D,  using  as  inputs  friction,  cavitation  and  lifting 
surface  theory.  The  friction  is  considerably  small  and  may 
usually  be  neglected,  and  the  pitch-diameter  ratio  may  be 
given  as, 

P/D  -  xxirx  (1  ♦  (A  (P/D) /  (P/D) )  *  tan(@i  +  c^  ))  (17) 
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where : 


1  +  (A  (P/D)) /(P/D)  =  (tan(Bi  +  o1  +  a2)0>7)/tan(8i  +  al) 

(18) 

*  57.3  Ka(0.0849  xC^  +  0.0152  *  a)  (in  degrees) 

a2  *  ab  +  (6i"  8)  *  (2/d  +  cos2Si(|-  D)  -  1) 

h  is  the  pitch  correction  factor,  and  is  the  vortex 
correction  factor  (may  be  neglected  due  to  its  magnitude) . 

is  a  function  of  the  cavitation  index  at  the  seventy 
percent  station,  oQ  ?,  which  is  in  turn, 

oQ  7  -  2  X  gx  a  X  J2/V*  (J2  +  4.84)  (19) 


where : 


H  is  the  absolute  pressure  at  the  shaft  centerline 
minus  the  cavity  pressure,  measured  in  feet  of  water. 

With  the  above  information  in  hand,  the  designer  can 
verify  the  shape  of  the  blade  by  determination  of  the  pro¬ 
peller  efficiency  and  the  delivered  power. 

A.  DESIGN  ITERATIONS 

In  the  previously  discussed  steps,  it  was  assumed  that 
the  ideal  thrust  coefficient,  CTi,  was  approximately  seven 
percent  higher  than  the  thrust  coefficient.  Equation  (9) . 


Once  the  initial  design  phase  has  been  concluded,  the  drag 
to  lift  ratio,  e,  can  be  determined  and  a  more  precise  value 
of  ideal  thrust  coefficient  may  be  calculated.  This  may 
be  used  for  additional  iterations  until  a  final  value  of 
thrust  coefficient  and  power  coefficient  have  been  obtained. 
The  ratio  of  thrust  coefficient  to  the  power  coefficient 
determines  the  efficiency  of  the  propeller. 

A  check  of  the  adsorbed  horsepower  is  made  to  see  if  the 
design  meets  the  requirements  of  power  available. 

SHP  (absorbed)  *  Cp  x  p  «  D2  «  x  S/(S50  x  8)  (20) 

B.  STRENGTH  OF  THE  BLADES 

The  final  design  step  is  the  determination  of  the  bending 
and  torsional  moments  of  the  blades. 


Figure  4.  Bending  and  Torsional  Moments 


From  the  geometry  of  Figures  4  and  5  and  from  structural 
considerations 
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“xO  1  “Tb  008  $i  +  MQb  3in  6i 


(21) 


M 


yG 

“Tb 


M, 


Qb 


-  Mrb  sin  Bi  "  MQb  cos  ei 


(22) 


1  3  2  f ^ 

JX  p  x  TT  x  RJ  x  v^/2  /  (x-xQ)  (1 -c  tan  Bi)-gji 

xo 


dx 


(23) 


r  x  p  X  it  X  R' 


x  V^/z  /  (tan  B)  (x-xQ) 


dCL, . 

x (1+e/tan  f^J-g—dx  (24) 


where  x  is  the  dimens ionless  station  being  analyzed. 

It  is  to  be  noted  from  Figures  4  and  5,  that  Mx0  and 
M^g  are  approximate  inasmuch  as  the  precise  expression  is 
a  function  of  4  rather  than  Bi»  The  difference  between  <j> 
and  B^^  is  ap,  where  aF  »  +  aj#  with  a  value  of  about  2  _ 

degrees.  Except  for  heavily  loaded  blades,  the  use  of  B^ 
instead  of  4>  is  valid,  since  with  heavy  loads  the  angle  of 
attack  criticality  increases. 

The  described  design  techniques  have  used  thrust  as  a 
base.  It  is  possible  to  use,  instead,  power  as  the  base, 
i.e.,  Cp  vice  C^. 
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Figure 


5.  Definition  Sketch 


III.  COMPUTER  AIDED  DESIGN 


A.  INTRODUCTION 

The  basis  of  a  Computer  Aided  Design  (CAD)  program  is 
to  make  available  to  the  designer  the  information  that  he 
needs  for  at  least  a  first  iteration.  To  be  efficient,  the 
CAD  should  have  the  following  characteristics: 

1.  Be  self-contained. 

The  entire  program  with  all  subroutines  should  be 
stored  in  a  single  data  source,  e.g.,  a  computer  disc,  so 
that  additional  programs  need  not  be  loaded  once  the  design 
process  has  commenced . 

2 .  Be  self-querying . 

The  variable  inputs  should  be  inserted  in  response  to 
a  computer  question  statement. 

3.  Be  labor  non- intensive. 

Fixed  parameters  should  be  included  in  the  program 
and  not  be  required  to  be  entered  each  time  the  program  is 
exercised. 

4.  Be  logical. 

The  computer  program  should  follow  a  logical  pro¬ 
gression  that  is  similar  to  the  steps  in  a  manual  calcula¬ 
tion,  unless  mathematical  considerations  dictate  otherwise. 

5.  Be  Segmented. 

The  program  should  contain  logical  breakpoints  to 


facilitate  future  editing,  including  changes,  deletions  and/ 
or  insertions. 


B.  APPROXIMATIONS 


The  use  of  empirical  data  requires  extensive  curve 
fitting,  interpolation  and  extrapolation.  Frequently  approxi¬ 
mations  may  be  made  to  the  data  that  signif icantly  reduces 
the  calculation  time.  Such  approximations  are  useful  and 
valid  only  if  they  provide  a  high  confidence  level  in  the 
results.  The  following  section  details  several  of  these 
approximations,  together  with  numerical  examples  as  a  proof 
of  confidence. 


An  expression  of  the  optimum  speed  of  advance. 


Jopt' 


was  developed  as. 


J1  "  (Cj  /J*  Jj-  0.56396  +  0.71765  x  Ln(Jx)) 


(Cj*-/J  +  0.71765/J1) 


where : 

J1  -  1.206/(cJ'5/J  +  0.63225) 

This  approximation  provides  a  confidence  level  of  98%. 
The  expected  efficiency  is  extracted  from  the  predicted 
performance  charts  extrapolated  from  40  knots  to  80  knots, 
for  blades  between  ten  and  fourteen  feet  in  diameter.  This 
extraction  was  performed  by  Lagrangian  curve  fitting.  To 
demonstrate  the  accuracy  of  these  approximations,  five  runs 
were  made  at  different  values  of  RPM  with  the  following 
conditions : 
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V  *  80  kts 


R  *  120,000  lbs 


z 


6 


SHP  =  45,000  Hp  w  *  t  *  0  p  »  2  slug/ft3 

Submergence  «  0.5 

The  results  of  these  runs  are  shown  in  Table  III  where 
the  subscript  'h'  refers  to  manual  calculations  and  the 
subscript  'c'  refers  to  computer  calculations.  All  calcula¬ 
tions  have  been  generated  by  this  author. 


TABLE  III 

SPEED  OF  ADVANCE/DIAMETER  (MANUAL  VS  COMPUTER) 


BUN 

BPM 

°h 

Jc 

yj 

opth 

jQptc 

*h 

nc 

Vc 

1 

400 

1.69 

1.69 

1.45 

1.46 

0.70 

0.70 

13.98 

13.91 

2 

500 

1.35 

1.35 

1.35 

1.36 

0.67 

0.70 

12.01 

11.92 

3 

600 

1.12 

1.13 

1.26 

1.28 

0.65 

0.65 

10.72 

10.56 

4 

700 

0.96 

0.97 

1.20 

1.21 

0.63 

0.62 

9.65 

9.57 

5 

800 

0.84 

0.84 

1.04 

1.15 

0.60 

0.59 

9.74 

8.82 

The  determination  of  the  uncorrected  radial  distribution 
of  the  blade's  characteristics  as  well  as  the  determination 
of  the  ideal  thrust  coefficient  requires  a  knowledge  of  the 
ideal  efficiency.  In  the  manual  mode,  values  were  extracted 
from  the  curves  of  Kramer's  Thrust  coefficient.  Once  again, 
curve  fitting  techniques  were  employed.  Because  of  the 


32 


shape  of  the  empirical  curves ,  the  numerical  approximation 
required  five  partitions. 


The  manual  and  computer  values  for  the  ideal  efficiency 
are  shown  in  Table  IV  for  five  runs  at  different  values  of 
rji  and  Cri  for  z  *  6. 


TABLE  IV 

IDEAL  EFFICIENCY  (MANUAL  VS  COMPUTER) 


Run 

X 

CTi 

nih 

nic 

1 

0.500 

0.0930 

0.941 

0.948 

2 

0.400 

0.1260 

0.930 

0.949 

3 

0.400 

0.1605 

0.927 

0.925 

4 

0.400 

0.1954 

0.905 

0^901 

5 

0.400 

0.2304 

0.875 

0.878 

The  camber 

correction 

,  coefficients 

,  k ^  and 

k^,  must  also 

be  approximated 

,  and  are 

based  on  the 

ratio  of 

expanded  blade 

area  to  disk  area,  A^/Aq .  Table  V  shows  a  comparison  of 
the  values  obtained  from  an  empirical  graph,  subscript  *g' 
to  the  computer  solution  values,  sucscript  'c'. 

The  graph  used  for  this  approximation,  [1] ,  has  a  maximum 
value  of  1.2  for  the  ratio  of  the  expanded  blade  area  to  disk 
area,  A^/A^.  Even  with  an  extension  of  this  value  to  nearly 
1.6,  the  computed  values  of  k^  and  k2  appear  to  be  well 
within  reason. 


TABUS  V 


CAMBER 

CORRECTION 

COEFFICIENTS 

(MANUAL  VS 

COMPUTER) 

*e/ao 

Klc 

K2g 

K2c 

1.54 

* 

1.06 

* 

2.31 

0.75 

1.03 

1.00 

2.10 

2.13 

0.72 

1.06 

1.06 

2.10 

2.09 

0.33 

1.01 

1.00 

1.70 

1.67 

0.25 

1.02 

1.00 

1.65 

1.67 

* 

Point  off  empirical  graph 

Two  additional  curve  fit  solutions  deal  with  the  cavitation 
index  and  the  lifting  surface  correction  angles  (a1  and  a2) . 

As  indicated  in  Equation  (16)  ,  the  correction  angles  are  a 
function  of  sigma  and  h.  To  indicate  the  degree  of  fit  of  the 
computer  generated  data  for  these  approximations.  Figures  6 
and  7  show  the  computer  plots  with  manually  extracted  data 
points  from  reference  6  superimposed. 

The  boldest  modification  in  the  development  of  this 
Computed  Aided  Design  program  was  the  complete  removal  of  the 
Goldstein  function  and  the  induction  factors  from  the  calcu¬ 
lations.  This  modification  is  required  in  order  to  simplify 
the  design  process.  The  Prandtl  solution  for  irrotational 
motion  of  a  screw  surface  in  an  inviscid,  incompressible, 
continuous  fluid  lends  itself  to  a  simplified  means  of 
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Figure  6.  Cavitation  Index  (Manual  vs  Computer)  [1] 
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Figure  7.  Lifting  Surface  Corrections  (Manual 
vs  Computer)  11] 


determining  the  circulation  and  the  lift  coefficient.  This 
approximate  solution  is. 


z  r  n 

IT  U  V 

a  a 


1  + 


r)  COS 


where : 

f  -  (1  +  Wq) 1/2(1  -  u/u0) 

p  *  Nx/Va 

«0  *  N/Va- 

However,  according  to  Goldstein  [5] , 

The  accuracy  of  the  (Prandtl)  approximation  increases 
with  the  number  of  blades  and  the  ratio  of  the  tip  speed 
to  the  velocity  of  advance,  but  for  given  values  of  these 
numbers,  we  have  found  no  means  of  estimating  the  error, 
since  the  exact  solution  has  not  been  found. 

Since  Goldstein's  work  in  1929,  NSRDC  [2]  has  developed 

a  more  precise  solution  utilizing  induction  factors. 

In  the  development  of  this  program  it  was  found  that  the 

direct  use  of  the  Prandtl  approximation  [5]  as  a  substitute 

for  the  more  complex  Goldstein  function  produced  excessive 

errors  when  compared  to  empirical  results.  It  was  therefore 

decided  to  use  a  modified  solution  which  did  not  contain  the 
2  2 

term  (u  /l  +  y  ) .  This  modified  term,  which  is  called  the 
Wilson  factor,  Wf,  provided  a  conservative  solution  with  much 
less  mathematical  manipulation  than  was  required  if  the 
Goldstein  function  of  Lerb's  induction  factor  method  were 
used.  The  Wilson  factor  is  defined  as. 
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(27) 


Wf  «.  (z/2ir)  x  cos”^  e"*^ 

Table  VI  shows  a  comparison  of  the  Wilson  factor,  Wf, 
and  the  Goldstein  function,  k,  at  several  radial  stations. 


Table  VI 

WILSON  FACTOR  VS  GOLDSTEIN  FUNCTION 


X 

K 

Wf 

ic/Wf 

0.400 

1.000 

1.209 

0.83 

0.475 

0.982 

1.160 

0.85 

0.550 

0.962 

1.101 

0.87 

0.625 

0.940 

1.031 

0.91 

0.700 

0.908 

0.946 

0.96 

0.775 

0.805 

0.841 

0.96 

0.850 

0.760 

0.704 

1.08 

0.925 

0.570 

0.511 

1.12 

PROGRAM 

DESCRIPTION 

1.  Segment  A 

The 

first 

segment  of 

the  program, 

lines  10  through 

630,  is  the  input  segment  where  the  desired  parameters  are 
inserted  in  response  to  queries  by  the  computer.  These 
inputs,  by  line  number,  are  shown  in  Table  VII. 
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TABLE  VII 


COMPUTER  PROGRAM  INPUTS 


Line 

Input  (units) 

190 

Ship  maximum  velocity  (kts) 

220 

Ship  maximum  resistance  (lbs) 

250 

Maximum  available  horsepower 

280 

Wake  fraction  (use  zero  if  V  * 

60  kts) 

310 

Thrust  deduction  factor 

340 

Number  of  blades 

370 

Fraction  of  submergence  of  propeller 

400 

Assumed  RPM 

430 

Fresh  or  salt  water  operation 

450 

Water  temperature  (Rankine) 

540 

Preliminary  propeller  diameter 

decision 

610 

Preliminary  propeller  diameter 

(ft) 

630 

3 

Fluid  density  (slugs/ft  ) 

2.  Sec 

jment  B 

The  second  segment  of  the  program,  lines  640  through 
930,  computes  and  prints  the  basic  parameters  for  the  ensuing 
calculations.  The  printout,  in  th».  computer  symbology, 
includes  the  thrust  load  coefficient  Ct,  from  existing 
formulae;  the  advance  coefficient,  J;  the  optimum  advance 
coefficient,  <11;  from  numerical  approximations;  the  optimum 
propeller  diameter,  Dl;  the  expected  efficiency,  E,  also  from 
numerical  approximations;  and  the  Shaft  Horsepower,  P. 


3.  Segment  C 


Because  of  the  complexity  and  non-linearity  of  the 
empirical  Kramer  Thrust  Coefficient  curves,  an  extensive 
curve  fitting  process  is  required  to  cover  the  entire  range. 
This  is  accomplished  in  a  multi-step  process  from  lines  940 
through  3170.  Lines  940  through  1370  are  used  to  approxi¬ 
mate  the  advance  ratio,  L,  in  order  to  enter  the  Kramer  curves, 
and  lines  1380  through  3170  determine  the  adjusted  value, 

Ln,  from  the  curves. 

4 .  Segment  D 

Once  the  adjusted  value  Ln  has  been  obtained,  the 
ideal  values  of  the  parameters  L^  Cfci  and  are  determined 
in  lines  3180  through  3300.  All  values,  except  are  from 
existing  equations. 

These  values  are  then  used  in  the  first  iteration  for 

the  ideal  thrust  coefficient,  Cti,  in  lines  3310  through  4460. 

» 

This  segment  is  concluded  with  a  printout  of  the 
delivered  ideal  thrust  coefficient,  ,  and  the  required 

ideal  thrust  coefficient,  for  comparison  purposes.  Also 

printed  at  this  point  are  the  values  of  the  tangent  of  the 
advance  angle,  tan  B,  the  values  of  the  Wilson  factor,  Wf, 
the  circulation,  G,  and  the  distribution  of  the  coefficient 
of  lift  times  the  chord-diameter  ratio.  All  calculations, 
with  the  exception  of  the  Wilson  factor,  are  from  existing 
equations. 
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5.  Segment  e 

Lines  4050  through  4460,  by  the  use  of  Simpson's 

Rule,  calculate  a  first  iteration  of  the  ideal  thrust  coeffi¬ 
cient,  and  a  final  iteration  of  the  blade  hydrodynamic 

delivered  ideal  thrust  coefficient,  Ct^,  as  well  as  the 
required  Cfc^. 

6 .  Segment  F 

A  second  iteration  of  the  ideal  thrust  coefficient  is 
accomplished  in  lines  4470  through  5020  in  order  to  verify 
the  blade's  hydrodynamic  pitch  angle.  This  is  necessary 
inasmuch  as  the  previous  calculations  were  based  on  the 
original  values  of  the  ideal  thrust  coefficient  vice  the 
iterated  values.  The  iteration  continues  until  the  delivered 
value  is  within  98%  of  the  required  value.  All  calculations 
are  based  on  existing  equations. 

7.  Segment  G 

Although  the  solution  may  be  in  hand  by  this  point, 
it  is  of  interest  to  recompute  lard>da  based  on  the  new  param¬ 
eter  values  and  then  to  use  this  value  as  an  alternate  method 
in  the  design  process.  This  is  available  in  lines  5040 
through  5390,  but  is  currently  by-passed  in  the  program. 

8.  Segment  H 

The  value  of  C^U/D)  is  recalculated  and  is 
printed  out  in  lines  5400  through  5430  with  the  aid  of  a 
subroutine,  and  the  shape  of  the  blade  is  determined  in  lines 
5440  through  5550,  in  addition  to  a  first  cut  approximation 
of  the  EAR  via  existing  formulae. 


rTi 


9 .  Segment  I 


The  geometry  of  the  foil,  i.e.,  cross-sectional  shape 
and  all  of  the  dimensions  at  each  predetermined  radial  sta¬ 
tion,  are  calculated  in  lines  5560  through  6300.  For  each 
station  there  is  a  printout  that  specifies  blade  element 
thickness ,  height  of  curvature  from  the  x-axis  and  blade 
width. 

10 .  Segment  J 

The  blade  pitch  is  corrected  for  cavitation,  friction 
and  the  .lifting  surface  effects  in  lines  6310  through  6570. 
The  printouts  for  this  segment  are  similar  to,  but  are  of 
corrected  values  of,  the  data  in  the  preceding  segment, 
with  the  final  pitch  diameter  ratio  being  determined. 

11.  Segment  k 

The  segment  from  lines  6600  through  6850  determines 
the  predicted  propeller  efficiency  and  manifests  that  the 
power  demanded  by  the  propeller  does  satisfy  the  available 
or  required  power. 

12.  Segment  L 

The  bending  and  torsional  moments  along  an  arbitrary 
radial  station  are  determined  in  the  segment  from  lines  6860 
through  7230  for  various  predetermined  stations. 

13.  Segment  m 

The  remainder  of  the  program,  through  line  8450,  is 
devoted  to  subroutines  that  perform  the  repetitive  calcula¬ 
tions  required  by  the  various  segments  and  data. 

41 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  program  is  an  approximate  method  for  aiding  in  the 
design  of  supercavitating  propeller  blades.  The  numerical 
approximations  complemented  by  the  conversion  factors  and 
other  constants  result  in  a  design  that  is  about  90%  accurate 
in  the  final  sense  as  compared  with  information  from  reference 
1.  Inasmuch  as  this  method  uses  essentially  the  same  formulae 
and  empirical  curves  as  the  'manual'  method,  one  could  not 
hope  for  much  better  accuracy. 

The  principal  advantage  to  this  CAD  program  is  that  it 
reduces  the  conceptual  design  phase  to  minutes.  It  also  pro¬ 
vides  a  straightforward  approach  to  the  problem  that  can  be: 

(a)  easily  followed,  and  (b)  modified  to  fit  the  specific 
needs  of  the  designer.  The  use  of  proven,  computerized 
equations  is  much  faster  than  table-look-ups  and  picking 
points  off  a  graph.  In  addition,  the  repeatability  of  output 
is  greatly  enhanced. 

The  program  will  definitively  demonstrate  what  combination 
of  fundamental  parameters  will  not  yield  a  satisfactory  design, 
e.g.,  negative  coefficient  of  lift  for  a  desired  forward 
motion.  This  information  provides  the  designer  with  great 
flexibility. 

If  the  program  does  produce  unsatisfactory  solutions,  the 
designer  may  vary  the  basic  parameters.  From  the  fluid 


dynamics  viewpoint  the  best  parameter  to  change  is  RPM, 
inasmuch  as  RPM  plays  a  major  role  in  the  design  of  propellers, 
and  by  adjusting  the  assumed  RPM  the  resulting  output  is 
changed . 

This  program  is  equipped  with  an  alternative  design  tech¬ 
nique  so  that  when  the  final  tangent,  8i#  at  x  *  0.7  is  veri¬ 
fied,  a  new  ideal  advance  ratio  may  be  calculated.  This 
alternative  is  not  a  portion  of  the  basic  program  but  has 
been  included  for  use  if  desired  at  a  later  time. 

A  comparison  with  a  manually  calculated  theoretical  design, 
as  shown  in  Tables  III,  IV  and  V,  demonstrates  that  this 
program  is  feasible  and  equally  as  accurate.  However,  the 
only  definitive  means  of  determining  the  validity  of  this 
program  is  to  make  a  comparison  with  an  experimental  design. 

This  program  can  be  used  for  all  supercavitating  propeller 
blades  in  an  extremely  wide  range  of  speeds,  submergences, 
thrusts  and  efficiencies. 

B .  RECOMMENDATIONS 

A  graphics  routine  would  enhance  the  program  and  might 
even  upgrade  this  program  to  a  Computer  Aided  Engineering 
(CAE)  and/or  Computer  Aided  Manufacturing  (CAM)  status. 

At  present  five*  and  seven-bladed  propellers  are  excluded 
in  this  program,  since  it  was  too  difficult  to  read  the  ad¬ 
vance  ratio  off  the  Kramer  Thrust  Curve.  This  could  be 
changed  in  the  future. 

The  effect  of  the  change  of  advance  ratio  on  the  final 
output  has  not  been  examined.  The  effect  of  this  alternative 
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should  be  challenged  in  the  future.  In  addition,  the  intro¬ 
duction  of  induction  factors  would  be  well  worth  the  effort. 
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280  INPUT  M 

290  PRINT  “THE  MAKE  FRACTION  IS“;W 

308  D ISP  ''Enter  thrust  deduction  factor  (see  progrto  explanation} 


310  INPUT  Td 

328  PRINT  -THE  THRUST  DEDUCTION  FACTOR  IS“;Td 
330  DISP  "Enttr  tht  nuwbvr  of  bltdtt  on  tht  pi 
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1371  GOTO  7120 

•380  FOR  X-.70  TO  .829  STEP  .879 
•398  GOTO  7130 
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APPENDIX  B 


SAMPLE  COMPUTER  OUTPUT 


The  pu rpose  of  this  prs^ru  is  to  4tttrtlM  tho  following  specification*  of  o 
•upercavitat ing  propeller: 

IlMttip  <os  a  function  of  on  assumed  PPM) 

*1  *0*  shop*  (Chord  dimensions  st  wsrious  rodiol  stktio ns/ 

Clad*  pitch  snd  pitch  distribution 
Foil  CkPbor  ond  thickness  distribution 
Vlad*  bonding  moment  *t  rodiol  station* 

In  order  to  give  you  the  desired  propeller  pleose  enter  the  known 
charocterist ics  of  the  ship  followed  by  pressing  CONT. 


Alpha  is  the  angle  of  attack 

Alpha!  is  the  correction  to  alpha  for  cavitation 

11  is  the  hydrodynamic  pitch  angle 

Cl  Is  the  coefficient  of  lift 

Cpc  is  the  final  pouer  coefficient 

Ct  is  the  thrust  load  coefficient 

Ctc  is  the  final  thrust  load  coefficient  - 

Cti  is  the  ideal  thrust  load  coefficient 

•  is  the  diameter  (ft./ 

d  is  the  percent  of  blade  width 

CAR  is  the  expanded  area  ratio  <A*xAo/ 

Ci  Is  the  ideal  propeller  efficiency 
Cpsiton  is  the  dragslift  ratio 
Cta  is  the  propeller's  efficiency 
C  is  the  circulation 

Kt  a  K2  are  the  camber  correction  coefficients 
L  ts  the  advance  ratio  (lambda/ 

Ln  I*  the  adjusted  advance  ratio 
LI  l*  the  Ideal  advance  ratio 
to  is  the  blade  width  at  »i,7 
I  ts  the  blade  width 
Amo  Is  the  moment  about  the  x  axis 
Aye  is  the  moment  about  the  y  axis 
t  ts  the  thickness  of  the  blade  element 
x  is  the  percent  radius  of  tn*  blade 
df  is  the  Hilton  factor 

y  is  the  height  of  curvature  of  the  blade  element 
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SNIP'S  MAX  SPECS  IS  SS 
SHIP'S  TOTAL  RESISTANCE  IS  IS0040 
SHIP'S  IMX  AVA1LA4LE  SMP  IS  4S0SS 
THE  HAKE  FRACTION  tS  • 

THE  THRUST  SEDUCTION  FACTOR  IS  S 
NUISER  OF  SLADES  OH  PROPELLER  IS  6 

THE  PERCENT  OR  FRACTION  OF  SUDMERCCNCE  EXPECTED  OF  THE  CRAFT  IS  .9 
ASSUMES  RPH  IS  <00 
CALCULATIONS  ARC  FOR  SU 

THE  ASSUMES  HATER  TEMPERATURE  IS  520.000400 
VISCOSITY  <HU>*  .MORI 3 


TN«  thrust  - 1  o«4  c*«ff  Icitnt  ■  .1474 

Th*  squat*  root  of  tM  thrutt'lttd  coefficient  is  .40*1 

Th*  «4«nk<  coefficient.  J  -  1.3912 

Tho  optloua  ttfvtMt  (MtflcttM  It  t.27*i 

Tho  optiowo  pr or*  1 1  or*  dltMttr  is  14.9434  ft. 

Tho  expected  efficiency  is  .4914 


Th«  SNP  for  this  ORttauo  prop»\1*r  is 

49244.402* 

L  •  .4M 

Cti*  .144 

Ln  -  .240 

Th*  idtal 

* 

A 

l 

*c 

<Ei>  ts 

.*29 

Li*  .433 

• 

X  - 

TAN  it 

C04  Si 

HP 

6 

C  I'D 

;4Ht 

I. 0019 

.474* 

1.2000 

.0049 

.02*4 

.4790 

.*107 

.73*3 

1.19*4 

.0074 

.0309 

•  9SM 

.7049 

.7040 

1.1012 

0 

.02** 

.4290 

.4*21 

.0223 

1.0312 

.4402 

.4241 

.7000 

.4144 

.4907 

.*443 

.4000 

.4294 

.7790 

.9942 

.0732 

.4407 

.0079 

.022Q 

.4900 

.3447 

.0912 

.7042 

.0044 

.0170 

.*290 

.4477 

.*090 

.3107 

.0049 

.0129 

THE  Sim  OF 

THE  PRODUCTS  IS  4.1 

IS? 

THE  SELIVCRES  Cti  • 

.122 

THE  REQUIRED 

Cti*. 140 

THE  VERIFYING  DELIVERED  Cti  • 

.2772 

THE  REQUIRED 

Cti  •  .1409 

THE  VERIFYING  DELIVERED  Cti  • 

.9940 

THE  REQUIRED 

Cti  -  .1409 

THE  FINAL 

TANGENT  OF 

DETA  SU4 

i  IS  .4204 

THAT  ANCLE  IN  DECREES  AT  0.7R  IS  II.R23S 

THE  FINAL  SETAt  HAS  BEEN  VERIFIES.  YOU  NAY  PROCEED  KITH  A  CONFIDENCE  FACTOR  OF 


Le«  1.6704 

THE  APPROXIMATE  (UITHIN  305!>  EXPANSES  MCA  RATIO  (EAR)  IS  .4940 


X 

1x10,7 

Ik 

Cl 

ALPHA 

ysiaax 

y*4K 

.400  l 

.0899 

1.8262 

.1690 

2.0000 

.0171 

.0313 

.479  1 

.0869 

1.0228 

.1766 

2.0000 

.0182 

.0331 

.990  1 

.0730 

1.8010 

.1793 

2.0080 

.0180 

.0324 

.(29  1 

.0460 

1.7597 

.1691 

2.0000 

.8170 

.0299 

.700  1 

.0000 

1.6704 

.1600 

2.0080 

.0156 

.0262 

.779 

.9030 

1.9156 

.  1534 

2.0000 

.0146 

.0221 

.090 

.7630 

1.2007 

.1472 

2.0000 

.0136 

.0174 

.929 

.9690 

.9403 

.1309 

2.0800 

.0123 

.0117 

O'l 

y 

»xl 

lift) 

.0079 

.0009 

.0023 

.0039 

.0129 

.0008 

.0033 

.0096 

.0900 

.0037 

.0087 

.0146 

.1000 

.0076 

.0138 

.0232 

.2000 

.0140 

.0221 

.0370 

.3000 

.0206 

.0200 

.0404 

.4000 

.0244 

.0344 

.0978 

.9000 

.0262 

.0399 

.0670 

.0000 

.0298 

.0497 

.0766 

.7000 

.0230 

.0922 

.0079 

.0000 

.0178 

.0984 

.0991 

.9000 

.0102 

.0640 

.1079 

.9900 

.0094 

.0664 

.1119 

i.oooo 

.0800 

.0606 

.1191 

EAR-  .3211 

Kl-  1.046  K2-  1.666 


Oxl 

y 

y'l 

<yxl>< 

y< 

.0079 

.0009 

.0002 

.0003 

.0008 

.0129 

.0008 

.0004 

.0806 

.0014 

.0500 

.0037 

.0016 

.0026 

.0062 

.1000 

.0076 

.0032 

.0093 

.0120 

.2000 

.0140 

.0062 

.0103 

.0249 

.3000 

.0206 

.0006 

.0143 

.0349 

.4000 

.0244 

.0102 

.0170 

.0418 

.9000 

.0262 

.0109 

.0192 

.0440 

.6000 

.0250 

.0107 

.0179 

.0432 

.7000 

.0230 

.0096 

.0160 

.0306 

.0000 

.0179 

.0074 

.0124 

.0299 

•  9000 

.0102 

.0042 

.0071 

.0171 

.9900 

.0054 

.0023 

.0039 

.0091 

1.0000 

.8000 

.0000 

.0000 

.0000 

X  TAN  81 

TAN  * 

3 

8i9.7 

•0.7 

C 

ALPHA 

ALPHA  1 

P'S 

.400  1 

.0019 

1.0179 

.7943 

.5939 

.5260 

.1698 

2.0008 

.0403 

1.4644 

.479 

.9107 

.6972 

.7006 

.5939 

.9260 

.1766 

2.0008 

.0408 

1.4693 

.990 

.7069 

.7403 

.6373 

.5839 

.9268 

.1793 

2.0000 

.0407 

1.4672 

.629 

.6921 

.6919 

.9774 

.5939 

.5260 

.1691 

2.0000 

.0402 

1.4697 

.700  ' 

.6100 

.9017 

.9260 

.9939 

.9268 

.1600 

2.0000 

.0395 

1.4727 

.779 

.9902 

.3294 

.4037 

.5539 

.9260 

.1534 

2.0000 

.0390 

1.4769 

.090 

.9009 

.4790 

.4467 

.9935 

.5268 

.1472 

2.0000 

.0309 

1.4817 

.929 

.4677 

.4402 

.4147 

.5939 

.5268 

.1399 

2.0000 

.0379 

1.40 61 

N  EPSILON 
.400  .070 

.470  .000 

.900  .000 

.CM  .070 
.700  .  070 
.770  .077 

.000  .077 

.000  .*70 
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LI-  .03 

THC  FOLLOWING  IS  RECALCULATES  USING  EPSILON  I 


X 

TAN  t< 

COS  SI 

HP 

C 

C  I'D 

.4000 

1.0819 

.6709 

1.2080 

.0069 

.8294 

.4790 

.9107 

.7393 

1.1996 

.0074 

.0303 

.9900 

.7069 

.7060 

1.1012 

.8000 

.0299 

.6290 

.6921 

.8223 

1.0312 

.0082 

.0281 

.7000 

.6180 

.8307 

.9463 

.8000 

.0234 

.7730 

.9902 

.0732 

.0407 

.0973 

.0220 

.0900 

.9089 

.9912 

.  7042 

.0066 

.0178 

.9290 

.4677 

.9098 

.9107 

.0049 

.0129 

THC  SUN  OF  THC  PRODUCTS  IS  4 
THC  DELIVERED  Cti  •  .122 

.807 

THC  REOUIRED 

Ctl  ■  .161 

THC  VERIFYING  DELIVERED  Ctl 

-  .2780 

THE  REOUIPED 

Cti  ■  .1688 

THC  VERIFYING  OELIVEPTD  Ct  I 

-  .9361 

THE  REOUIRED 

Cti  >  .1608 

THE  FINAL  TANGENT  OF  SETA  SUS  i  IS  .6206 
THAT  ANGLE  IN  DEGREES  AT  0.7R  IS  31.S2S0 

THC  FINAL  BETA!  HAS  SEEN  VERIFIES.  YOU  NAY  PROCEED  UITH  A  CONFIDENCE  FACTOR  OF 

SIX. 

etc-  . 1140  C#<-  .14*4  ETA-  . 7923 

SHP<ae>orb*d>«  212. 9914  AVAILABLE  SHP  IS  49000.0000  (ih*  propdltr  cin  net 
deeawd  sort  than  the  aval  tab) •  SNA) 

THC  FOLLOWING  CROUP  IS  FOR  Xe-.4O0  (NOTE  ALL  NONENTS  ARC  IN  FT  LBS> 

AT  Re-  .4000  THC  SENDING  nonCNT  IS  S0S3S.S6S7 
AT  Xe-  .4000  THE  TORSIONAL  NONCNT  IS  34113. 8SI0 
X  Hxo  Nyo 

.401  99407. S102  1309S.SS29 

.479  60300.0900  0793.9003 

.991  ■  60894.3779  4373.0044 

.029  61007.9100  692.8379 

.'700  60999.9444  -2912.9919 

.779  60707.7914  -9218.0882 

.090  60943.0929  -7344.0786 

.929  60298.1747  -9396.7943 

THC  FOLLOWING  GROUP  IS  FOR  No-  .4790 
AT  Xe-  .4730  THC  BENDING  NONCNT  IS  83811.3339 
AT  Xe-  .4730  THC  TORSIONAL  NONCNT  IS  99849.1160 
X  Axe  Nyo 

.479  103421.6972  16141.8943' 

.990  106311.7099  .  8491.9793 

.623  106631.7189  1990.7739 

.700  106991.4329  -3943.0232 

.779  106328.8192  -8274.3030 

.030  109933.3317  -12343.7920 

.929  109463.9400  -13863.8093 

THC  FOLLOWING  GROUP  IS  FOR  Xe*  .9900 
AT  Xe-  .9300  THC  SCHDIHG  NONCNT  IS  114027.9004 
AT  Xe-  .9900  THC  TORSIONAL  NONCNT  IS  74399.7940 


X 

Nxe 

Hye 

.990 

136371.9960 

12399.4846 

.620 

136071 . 7231 

4013.8192 

.700 

136099.6694 

-3091.2930 

.779 

136623.2034 

•9169.4928 

.000 

136171.3779 

-14399.1404 

.929 

133616.1224 

•18927.4090 

THC  FOLLOWING  CROUP 

IS  FOR  Xe-  . 

6230 

AT  Xe- 

.6290  THE  SENDING  NONCNT 

IS  128984.6674 

AT  Re¬ 

.6290  THC  TORSIONAL  NONCNT  IS  80823.8079 

ft 

Nxe 

Aye 

.020 

191727.0247 

6721.9393 

.700 

191872.2420 

•1197.4030 

.770 

191670. 9019 

•7902.8091 

.000 

191296.9700 

•13710.1237 

.920 

190719.0903 

•10741.3019 
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TMt  FOLLOWING  GROUP  IS  FOR  Xo*  .7000 
AT  Xo*  .7000  THE  SENSING  NONENT  IS  130081.8343 
AT  Xo*  .7000  THE  TORSIONAL  MOHENT  IS  77733.3941 
X  Nxo  Nyo 

.700  131521. 1627  2250. 6460 

.775  151471.9207  -4474.3349 

.050  151189.1292  -10270.7420 

.925  1 387*2. 8676  -13307.3707 

THE  FOLLOWING  GROUP  IS  FOR  Xo*  .7730 
AT  Xo*  .7730  THE  SENSING  MOMENT  IS  119319.4017 
AT  Xo*  .7750  THE  TORSIONAL  MOMENT  IS  63324.3146 


X 

Hxo 

Hyo 

.773 

134026.2339 

1113.3300 

.830 

133969.0351 

-4090.9931 

.923 

133737.0394 

-0626.2366 

THE  FOLLOWING  GROUP  IS  FOR  Xo*  .0300 
AT  Xo*  .0500  THE  SENSING  MOMENT  IS  96297.3690 
AT  Xo-  .0300  THE  TORSIONAL  MOMENT  IS  43396.6491 
X  Nxo  Hyo 

.030  103396.34 76  4023.1160 

.923  105690.4106  1302.7247 

THE  FOLLOWING  GROUP  IS  FOR  Xo*  .9230 
AT  Xo-  .9230  THE  SENSING  MOHENT  IS  61013.7364 
NT  Xo*  .9230  THE  TORSIONAL  MOMENT  IS  12330.3974 
X  Nxo  Nyo 

.923  60306.9012  14479.3034 

THE  FOLLOWING  GROUP  IS  FOR  Xo*  1.0000. 

THERE  IS  NO  GROUP  FOR  Xo*1.0l  THIS  RUN  IS  COMPLETE. 
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APPENDIX  C 


SAMPLE  CALCULATIONS 


Calculations  have  been  made  at  the  seventy  percent  radial 
station  (x  *  0.7)  in  both  manual  and  computer  modes  with  the 
following  input  data: 


Maximum  Velocity 

80  kts  {Design  Point) 

Power  Available 

45,000  SHP 

Ship  Resistance 

120,000  lbs 

Propeller  Speed 

600  RPM 

Submergence 

0.5 

Wake  Fraction 

0.0 

Number  of  Blades 

6 

Density 

2.0  (Salt  Water) 

.  CALCULATIONS 

LINE  PARAMETER 

COMPUTER 

MANUAL 

570 

D 

9.17 

9.2  ft 

640 

Va 

80  kts 

80  kts 

650 

T 

120,000  lbs 

120,000  lbs 

660 

Ct 

0.1162 

0.1162 

720 

J 

1.126 

1.126 

740 

Jopt  “  J1 

1.279 

1.28 

790 

D1 

10.563  ft 

10.56  ft 

880 

L 

0.400 

0.407 

890 

Cto 

0.150 

0.153 

900 

Cti 

0.160 

0.164 
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LINE 

PARAMETER 

COMPUTER 

MANUAL 

3250 

Ei 

0.975 

0.922 

3290 

Li 

0.433 

0.4414 

3420 

Tbi 

0.618 

0.6306 

3430 

Tb 

0.5817 

0.582 

3440 

B 

0.527 

0.527 

3720 

Mou 

0.0518 

0.0518 

3730 

Mouo 

0.00740 

0.00740 

3770 

Wf 

0.9463 

0.9463 

4350 

Ctidl 

0.122 

0.1675 

Note:  The  manual  value  is 
computer  program  is  based 

based  on  but  one 
on  ten  points. 

point,  while 

5360 

Circ 

0.0080 

0.0084 

5370 

Clldl 

0.0254 

0.0266 

5420 

L0 

1.678  ft 

1.756  ft 

5480 

C1 

0.16 

0.16 

5710 

Y 

lmax 

0.0156 

0.0156 

5720 

Y 

max 

0.0230 

0.0274 

5820 

Y 

0.0230 

0.0274 

5830 

Tlr 

0.0522 

0.0399 

5960 

Ear 

0.3211 

0.3200 

6150 

K2 

1.6666 

1.70 

6160 

K1 

1.000 

1.04 

6240 

Y1 

0.0096 

0.0109 

6250 

Ylc 

0.0160 

0.0193 

6260 

Yc 

0.0386 

0.0339 
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LINE 

PARAMETER 

COMPUTER 

MANUAL 

6500 

Alpha ^ 

0.0395 

0.0394 

6560 

Pdf 

1.4727 

1.56 

6690 

Eps7 

0.078 

0.078 

6720 

Ctc 

0.116 

0.155 

6820 

Eta 

0.79 

0.78 
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